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Abstract 

We investigate the process pp —>■ HZ^ + A at the -/s = 14 TeV LHC up to the QCD next- 
to-leading order (NLO), and discuss the kinematic distributions of final products after on-shell 
Higgs and Z-boson decays by adopting the narrow width approximation. The dependence of 
the leading order (LO) and the QCD NLO corrected integrated cross sections on the factoriza¬ 
tion/renormalization scale is studied. Our results show that the LO integrated cross section and 
kinematic distributions are significantly enhanced by the NLO QCD corrections, and the NLO 
QCD A-factor strongly depends on the observables and phase space. We conclude that in preci¬ 
sion experimental data analyse for probing the HZ^ coupling we should consider the NLO QCD 
corrections and put proper constraints on lepton-pair invariant mass to reduce the background. 

PACS: 11.15.-q, 12.15.-y, 12.38.Bx 
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1 Introduction 


In the standard model (SM) the Higgs boson serves for the breaking of the electroweak symmetry 
and the generation of the fundamental particle masses mm- Studying the Higgs mechanism is one 
of the main goals of the CERN Large Hadron Collider (LHC). In 2012 both the ATLAS and CMS 
collaborations announced the discovery of a new boson, whose properties are relatively close to the 
long awaited SM Higgs boson with mass of mn ~ 126 GeV [311]. After the discovery of the Higgs 
boson, the main task of further experiments is to determine its properties. Particularly, the precise 
determination of the Higgs boson couplings is imperative for verifying the validity of the SM and the 
existence of new physics at high energy scale. 

The LHC first runs at 7 and 8 TeV are completed, and the latest data from them show that all 
the properties ISHS] of the new boson measured so far are well consistent with that of the SM Higgs 
boson, but it is well known that there are some theoretical difficulties associated with the SM Higgs 
sector. For example, the famous hierarchy problem, which is associated with the quadratic radiative 
corrections to the SM Higgs mass, is one of the difficulties, and there is no way to solve this problem 
in the SM. So new physics effects are still expected to solve these difficulties. There are many model 
candidates of new physics predicting sizeable deviations from the the Higgs couplings in the SM. In 
the SM some Higgs couplings of the types guYV and gnvv'V" are absent at the tree-level, and they 
would be particularly sensitive to new physics mm- With the increases of the LHC luminosity 
and colliding energy, we can collect statistically enough events for most of the important multi-body 
production processes. Obviously, precision measurements require accurate theoretical predictions for 
both signal and background. In the last few years, the phenomenological results including the next-to- 
leading order (NLO) QCD corrections to the Higgs boson production associated with di-gauge-boson 
at the LHC, such as pp —^ HWW, pp —)• HW^-y, pp —)• HW^Z, have been studied [HHHj- 

The HZ^ production at the LHC also offers the possibility to directly investigate the HZ-j, HZZ-j 
and HZ'j'y anomalous Higgs gauge couplings [ISlilS]) as they would cause deviations from the SM 
predictions. Moreover, an accurate estimate of the pp —)• HZ^ + X process followed with subsequent 
Higgs and Z-boson decays could provide the direct observable predictions in searching for possible new 


2 


physics. With the transverse momentum and rapidity cuts of = 5 ~ 20 GeV and = 2.5 on 

the final photon, the cross section for the HZ’j production at the 14 TeV LHC is about 2.4 G.O fb, 
and therefore only about 170 ~ 420 fb~^ integrated luminosity is required to produce about 1000 
events. 

In this paper, we make a precision calculation for the pp —)• HZ^ + X process at the LHC including 
the NLO QCD corrections with on-shell Higgs and Z-boson decays in the narrow width approximation 
(NWA), but do not provide the detailed strategy to extract the information of the anomalous Higgs 
gauge couplings from the pp —)• HZ^+X process. In section H we give the description of the analytical 
calculations for the LO cross section and the NLO QCD radiative corrections to the pp —)• HZ^ + X 
process. In section HI we present some numerical results and discussions. Finally, a short summary 
is given. 

2 Description of the computation 

In this section we describe the analytical calculations at the LO and QCD NLO for the pp —)■ HZ^ + X 
process. 

A. LO calculation 

In our calculations we employ FeynArts 3.4 package m to generate LO and QCD NLO Feynman 
diagrams and their corresponding amplitudes. The algebraic manipulations on the amplitudes are 
implemented by applying FormCalc 5.4 programs [TB]. We neglect the masses of u-, d-, c-, s-quarks. 
Due to the smalless of (anti)bottom-quark density in proton, the LO contribution to the cross section 
from the pp —?■ 66 — HZ'j + X process at the 14 TeV LHC is less than 0.6%. Therefore, we do not 
consider the partonic process of 66 annihilation in our calculations. Then the contributions to the 
cross section for the parent process pp — HZ^ + X come from the following partonic processes, 

qipi) + q{P 2 ) ^ H{p 3) + Z{P4)+-f{p5), {q = u,d,c,s), (2.1) 

where pi, p 2 and ps, p 4 , ps represent the four-momenta of the incoming partons and the outgoing 
H, Z and photon, respectively. The LO Feynman diagrams for the partonic process qq —?• HZ^ are 
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Figure 1: The LO Feynman diagrams for the partonic process qq —)• HZ^. 


shown in Fig(TJ Since the Yukawa coupling strength is proportional to fermion mass and the u-, d-, c-, 
s-quarks are considered as massless, there is no contribution from the Feynman diagrams with internal 
Higgs boson line and the Higgs emission from initial quark. The LO amplitude for the qq —)• HZ^ 
partonic process involves QED soft and collinear IR singualrities since the photon is radiated from 
massless quark. To avoid these QED IR singularities and obtain an IR-safe LO result, we take the 
transverse momentum and rapidity cuts on the final photon [Prcut^ IVcutl) declared in the following 
section. The LO matrix element for the partonic process qq —)• HZj can be expressed as 


■Mlo — -Mto + -^LO 

o3 


m Qqtnz_y ^ n/ q q b\ (^1 i' { \ ^ 


+ 


ie^Qqmz_, , y{h-h) ur q 9 5^ ^ ^ 
-v{p2h — — y9v - 9 a1 )u{pi) 


c2 ^2 


{pb -P 2 Y 




(2.2) 


where = sin 9w, 0 ^, = cos 6\y, 9y = ^Tg—Qq sin^ dw, = \Tg, Tg and Qg are the third component 
of weak isospin and the electric charge of quark q separately. 

The LO cross section for the partonic process qq —)• HZ^ can be obtained by performing the 
integration over the phase space expressed as below, 

where p is the three-momentum of one initial parton in the center-of-mass system (c.m.s), y/S is the 
colliding energy in partonic c.m.s, the summation is taken over the spins and colors of the initial 
and final states, and the bar over the summation indicates the averaging over the intrinsic degrees of 
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freedom of initial partons. is the three-body phase space element defined as 


dQs = ipi+P2-'^Pi] n 


d^Pi 


i=3 / j=3 


(27r)32& 


(2.4) 


By convoluting with the parton distribution functions (PDFs) of the colliding protons, we 


obtain the LO total cross section for the parent process pp —)• HZ’j -|- X as 

aLO = 'y] / dxidX2 Gq/p^{xi,Pf)Gg/p^{x2,Pf)dqjV§ = XiX2'/s) + {1 ^ 2) 

„ Jo I 


(2.5) 


where Gg/p represents the PDF of parton q in proton P, xt {i = 1, 2) describes the momentum fraction 
of a parton in proton, y/s is the colliding energy in the rest frame of proton-proton system, and p/ is 
the factorization scale. 


B. NLO calculation 

In the NLO calculations we use the dimensional regularization (DR) method in D = 4 — 2e 
dimensions to regularize the ultraviolet (UV) and infrared (IR) divergences. The NLO QCD corrections 
to the pp —7- HZ^ + X process are constituted distinctly by the following three parts: (1) the virtual 
correction, (2) the real gluon and light-(anti)quark emission corrections, (3) the collinear counterterms 
of the PDFs. The virtual NLO QCD correction to the qq —)• HZ^ partonic process consists of self¬ 
energy, vertex, box and counterterm diagrams. The one-loop Feynman diagrams are shown in Figl2j 
We follow the definitions of tensor and scalar one-loop integral functions in Refs. |19ll20j . and use 
the Passarino-Veltman (PV) method |19ll21j to reduce tensor integrals to the linear combinations of 
tensor structures and coefficients, where the tensor structures depend on the external momenta and 
the metric tensor, while the coefficients depend on scalar integrals and kinematic invariants. The whole 
reduction manipulations of a tensor integral to the lower-rank tensors and further to scalar integrals, 
is done and numerically calculated by using the LoopTools-2.2 library [18] and the FF package |22j . 

In the virtual correction calculation we need the wave function renormalization constants for quark 
fields. We introduce the renormalization constants SZ^^ ^ ^ for massless quark {q = u, d, c, s) fields 
dehned as 

'^q,L,R = (1 + (2-6) 
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Figure 2: The one-loop Feynman diagrams for the partonic process qq —?■ HZ'j. 


In the modified minimal subtraction (MS) renormalization scheme the renormalization constants for 
the massless quarks are expressed as 

Q/ 

- ^ir), (2.7) 

where Cp = 4/3, Apv = - 7 e + ln(47r) and Ajr = -'jp + ln(47r). 

After the reduction for tensor integrals, the amplitude for loop corrections involving one-loop scalar 
integrals contains both UV and IR divergences. The UV divergence is vanished after performing the 
renormalization procedure. But the total QCD NLO amplitude for the subprocess qq — HZj still 
contains QCD soft/collinear IR singularities. We adopt the expressions in Ref. |23j to deal with the 
QCD IR divergences in Feynman integrals, and apply the expressions in Refs. |24H26j to implement 
the numerical evaluations for the QCD IR-finite parts of Wpoint scalar integrals. According to the 
Kinoshita-Lee-Nauenberg (KLN) theorem |27j . these IR singularities will be cancelled by adding the 
contributions of the real gluon/light-(anti)quark emission subprocesses, and redehning the PDFs at 
the QCD NLO. 

Since we put the transverse momentum, rapidity cuts on the final photon and a resolution cut on 
the photon and final jet throughout our LO and QCD NLO calculations, the numerical results for 
the real gluon/light-(anti)quark emission subprocesses are QED IR safe. While the real gluon/light- 
quark emission processes contain the QCD soft and collinear IR singularities. Technically, we isolate 
the QCD soft and collinear IR singularities by adopting the two cutoff phase space slicing (TCPSS) 
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Figure 3: The Feynman diagrams for the real gluon emission subprocess qq —)• HZ'j + g. 



Figure 4: The Feynman diagrams for the real light-quark emission subprocess qg —)• HZ^ + q. 

method [28| . The Feynman diagrams for the real gluon and light-quark emission subprocesses are 
shown in FiglH] and FigHJ separately. Before our numerical calculations, we checked the UV and IR 
divergence cancelations both analytically and numerically. To verify the implementation of the TCPSS 
in right way, the independence of the NLO QCD corrected total cross section on the soft cutoff 5s 
are checked in the range of 1 x 10“® < 5^ < 1 x 10“^ with 6c = (5s/50. In the further numerical 
calculations, we fix = 1 x 10“^ and 5c = 2 x 10“^. Furthermore, we numerically compared our 
NLO QCD corrected cross sections with those obtained by using program MadGraph5_aMC@NLO of 
version 2.2.2 [5U], and find they are in good agreement with each other within the Monte Carlo errors. 
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3 Results and discussions 


In this section we present and discuss the numerical results for the HZ^ associated production at 

the Vi = 14 TeV LHC at both the LO and QCD NLO. We use the CTEQ6L1 and CTlOnlo PDFs 

in the LO and NLO calculations, respectively. The strong coupling constant is determined by taking 

one-loop and two-loop running as(/x) for the LO and NLO calculations separately, and setting the 

QCD parameter as Nf = 5, = 165 MeV for the CTEQ6L1 and = 226 MeV for the 

CTlOnlo. We set the factorization and renormalization scales to be equal, and take /r = /rj = /r,, = /^o 

by default unless stated otherwise. The central scale is defined as fxo = Et/2 = where 

i 

Ex^i = and the summation is taken over all the transverse energies of final particles. The 

related SM input parameters are taken as [31I41I30] 

a-^ = 137.036, mw = 80.385 GeV, mz = 91.1876 GeV, mn = 126 GeV. (3.1) 

To strip the QED soft and collinear IR singularities at both the LO and QCD NLO, we put the 
following transverse momentum, rapidity and resolution cuts on the final photon, i.e.. 


Pt > Pt,cuV 

1 ^ 7 ! 

< \y2ut\ = 2 . 5 , 

(3.2) 

fi'yj (5q 

or 

j O' ^ COS 

P^T^Pt 1 r • 

1 — COS Oo 

(3.3) 


where 6q is a fixed separation parameter which is set to be 0.7. The condition of Eq. (l3.3h implies that 
the final jet can arbitrarily close to the photon as long as the jet is soft enough. In this way, we can 
preserve the full QCD singularities, which cancels against the virtual part, but it does not introduce 
divergence from the interaction between photon and massless quark-jet m- The limitation in Eq. (13.31) 
is to remove the QED collinear IR singularity due to a photon radiated from a final light-quark-jet j 
in the NLO calculation for the real light-quark emission processes. Then we accept the ElZ^ + jet 
events only if all the limitations in (13.2p and (13.3p are satisfied. 

In Figl5] we display the renormalization/factorization scale dependence of the LO, NLO QCD 
corrected total cross sections and the corresponding iL-factor for the pp — )• HZ^ + X process at the 
y/s = 14 TeV LHC by setting pr = fJ-f = tJ'- The LO and NLO QCD corrected integrated cross 
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Figure 5: The dependence of the LO, NLO QCD corrected total cross sections and the corresponding 
iL-factor for the pp —>• HZ^ + X process on the factorization/renormalization scale. Here we assume 
p = Pr = fJ'f and define the central scale as po = Et/2. 


PLut (GeV) 

O-LO (/&) 

ctnlo ifb) 

K 

5 

4.780(4) 

5.99(3) 

1.25 

10 

3.105(3) 

4.01(2) 

1.29 

15 

2.287(2) 

3.00(1) 

1.31 

20 

1.785(1) 

2.37(1) 

1.33 


Table 1: The LO, NLO QCD corrected integrated cross sections and the corresponding iL-factors with 
different cuts on p^ for the HZ^ production at the 14 TeV LHC. 


sections are 1.785 fb and 2.37 fb, respectively, and the corresponding iL-factor is 1.33 at the central 
scale p = pq. If we defined the relative scale uncertainty as rj = [max{a{p)) — min{a{p))] /cj(^o) with 
p e [0.25/io, 4:/io], we get rj = 9.85% and 4.24% for the LO and NLO QCD corrected corss sections, 
respectively. We can see that the dependence of the NLO QCD corrected total cross section on the 
factorization/renormalization scale is significantly reduced compared with that of the LO integrated 
cross section. This makes the theoretical predictions much more reliable. 

We present the LO, NLO QCD corrected integrated cross sections and the corresponding iL-factors 
for different cuts on p^ in Tab{TJ We can see that the LO and NLO QCD corrected total cross sections 
strongly depend on while the iL-factor behaves not so sensitively to the cut. As shown in TablH 

we can get a sizeable decrease in the total cross section with the increase of p^ for the pp —)• HZ^+X 
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p' (GeV) 


Figure 6: The LO, NLO QCD corrected transverse momentum distributions of the final particles and 
the corresponding i^-factors for the pp —)■ HZj + X process at the LHC. (a) distributions, (b) p^ 
distributions, (c) pj distributions. 


process. In following numerical calculations we fix = 20 GeV as the default choice. 

We depict the LO, NLO QCD corrected distributions of the transverse momenta and the corre¬ 
sponding iL-factors for the Higgs, Z-boson and photon produced by the pp —)■ HZ^-\-X process at the 
y/s = 14 TeV LHC in FigslU^a), (b) and (c), respectively. In FigsEJ^a) and (b), the curves for and 
at the LO and QCD NLO peak at the position of pT ~ 50 GeV, and their iL-factors are 1.31 and 
1.33, respectively. FigEJ^c) shows that both the LO and NLO QGD corrected p'^ distributions decrease 
rapidly with the increase of the transverse momentum of photon. We can see from FigsE^a,b,c) that 
the transverse momentum distributions for the Higgs, Z-boson and photon {daio/dp^, daio/dp^, 
daio/dp'^) are significantly enhanced by the NLO QCD corrections. 

The final photon can be directly detected in experiment, while the produced on-shell Higgs and 
Z-boson are unstable particles and can be detected by their decay products. In order to investigate 
the kinematic distributions of final directly detected particles, we apply the NWA in analysing the 


10 
















































differential cross sections of the final Higgs and Z-boson decay products. We choose the Higgs boson 
decay channel of H —)• r+T“ with nir = 1776.82 MeV [3U] and Z-boson decay channel of Z —>• 

{i = e, fi) as the Higgs and Z-boson signals separately. By adopting HDECAY program |32] with 
input parameters from Ref. [3U], we get Br{H —)• r+r“) = 5.897%, and take Br{Z —>■ = 

Br{Z e+e") + Br{Z —)• = 3.363% + 3.366% = 6.729% [30]. Then the signature for the 

HZy production including the subsequent decays at the LHC can be written as 

pp ^ HZj ^ + X (£ = e,/i). (3.4) 

This signal is detected as an event including one r-pair, one e(/r)-pair and a photon. For photon 
separation with other final particles in the rapidity-azimuthal angle plane, we impose the R cuts 
between photon and other final particles as below: 

Rr'-f > 0.4, Ri^ > 0.4, Rj^ > 0.7, (3.5) 

where i = e, p, and j denotes a jet with transverse momentum plp > 30 GeV. 

The final photon in the signal process pp —)• HZy —)• {I = e, p) is only emitted 

from initial parton. However, the HZ associated production followed by the subsequent decays of 
H —)• T+r“ 7 , Z —>■ or 77 —)• r+r", Z also leads to the event and contributes 

at the same order as the signal process. This process, denoted as 

pp ^ HZ ^ + X {i = e,p), (3.6) 

is the main background in measuring the HZy coupling via the pp — HZj —?■ + X 

process. We estimate the background process pp ^ HZ ^ at the LO in the NWA by 

adopting CTEQ6L1 PDF, and take Tz = 2.4952 GeV [30]. The total decay width of SM Higgs boson 
is obtained by using HDECAY program as Th = 4.38 x 10“^ GeV. In analysing the 
events, we adopt the event selection criteria shown in Eas. (|3.2|) and (|3.5I) . and impose the following 
invariant mass constraints on the final lepton pairs to suppress the background contribution: 
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< mu + A 


niH — A < M^+^- 

mz - ^< Mi+i- <mz + A, {I = e, /i), 

Then we obtain the background over signal as 

(^Lo{pp —^ HZ —^ T+r“^+£“7) _ r 8.1%, 

O'LoiPP HZ'f —)> T+T~l+i~'^) \ 1.5%, 

Therefore, we can conclude that the background events with photon radiated from final charged leptons 

can be reduced distinctly in probing the HZj coupling by taking proper invariant mass constraints 

on T-pair and e(/z)-pair. In the following, we neglect the background contribution and only consider 

the signal process pp —)■ HZj —)■ in investigating the kinematic distributions of the final 

produced leptons. 

In FigslTl^a) and (c) we present the LO, NLO QCD corrected transverse momentum distributions of 
and positively charged lepton (i = e,/r), and the corresponding iiT-factors at the ^/s = 14 TeV 
LHC, respectively. As shown in Figsd^a) and (c), the QCD corrections always enhance the LO 
differential cross sections daio / dp'^ and daLO / dp^ ■ Both the LO and NLO QCD corrected transverse 
momentum distributions of r'*' and reach their maxima at the positions of pTj^ ~ 40 GeV with 
K = 1.31 and ~ 30 GeV with K = 1.31, respectively. Figsd^b) and (d) are for the LO and 
NLO QGD corrected rapidity distributions of r"*" and f"*" separately. Both yH and reach their 
maxima at y = 0, with their AT-factors being around 1.36. We can see from all these four figures that 
the NLO QGD corrections do not change the line-shapes of the transverse momentum and rapidity 
distributions, while enhance the LO differential cross sections significantly in all the plotted kinematic 
regions. With the transverse momentum and rapidity cuts of plp ,cut = 5 ~ 20 GeV and y^^^ = 2.5 
on the final photon, the accumulated luminosity of 1000 ~ 2660 fb~^ is required to produce 25 
events via pp —)• HZj —+ X channel at the 14 TeV LHG. 

4 Summary 

In this paper we investigate the NLO QGD corrections to the HZ'j production followed by subsequent 
Higgs and Z-boson decays at the ^/s = 14 TeV LHC. We study the dependence of the LO and NLO 


(A = 5 or 10 GeV). (3.7) 


A = 10 GeV 
A = 5 GeV ■ 


(3.8) 
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Figure 7: The LO, NLO QCD corrected transverse momentum and rapidity distributions of r'*' and 
the positively charged lepton, and corresponding iil-factors for the pp —)• HZ^ + X —)• 'yX 

processes at the y/s = 14 TeV LHC. (a) pTj^ distributions, (b) distributions, (c) p^ {i = e, p) 
distributions, (d) {i = e, p) distributions. 
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QCD corrected cross sections on the factorization/renormalization scale, and our resnlts show that 
the scale uncertainty of the NLO QCD corrected cross section is reduced compared with that of the 
LO cross section. We present the LO and NLO QCD corrected distributions of transverse momenta 
and rapidities of the decay products of Higgs and Z-boson. We find that the NLO QCD radiative 
corrections are significant, and notably modify the LO kinematic distributions. We see also that the 
iL-factor is distinctly related to phase space region and kinematic observable. We conclude that the 
NLO QCD corrections should be considered in precision experimental data analyse in measuring the 
pp —)• HZ^ + X process, and the background events with a photon radiated from final charged lepton 
can be reduced in probing the HZ^ coupling by putting proper invariant mass constraints on final 
r-pair and e(//)-pair. 
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